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A correlation analysis of a thickness tendency equation including 
heating and vertical motion terms was carried out. One of the main pur- 
poses of the investigation is to provide a firmer foundation for a 1000- 
500 mb thickness forecast scheme. 

A large data sample was used and single and multiple linear coeffi- 
cients were calculated. These were computed for a wintertime situation 
with separate analyses for Atlantic and Pacific areas. 

The resulting correlation coefficients for the advection, vertical 
motion and heating terms gave substance to the hypothesis that the observed 
change in thickness may be fairly well approximated by the change due to 
horizontal advection only. The vertical motion term gave a smaller but 
significant correlation, while heat exchange terms did not correlate 
significantly. 

The author wishes to express his appreciation to Professor G. J. 
Haltiner for his guidance and encouragement, to Mrs. Jean Bow of the 
computer center at the U. S. Naval Postgraduate School, for her help in 
Fortran programming, and especially to Mr. James R. Clark of Meteorology 
International, Inc., for his assistance in data collection and map re- 
production. The author also wishes to express appreciation to the U. S. 
Naval Fleet Numerical Weather Facility for making their data and computer 


facility available for this study. 
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i. INTRODUCTION 

The ability to produce analyses and prognoses of weather maps by 
numerical methods using high speed digital computers on an operational 
basis is now a reality. Most of the efforts in numerical weather pre- 
diction have been in relation to the pressure and temperature fields; 
however some pioneering attempts have been made to predict precipitation, 
clouds and visibility. 

Perhaps the most successful forecasts have been of the 500-mb 
contours. However, somewhat less success has been achieved in obtaining 
a good surface or 1000-mb prognosis. Many schemes have been tried with 
varying degrees of success. One method that has a great deal of promise 
makes use of the 500-mb and 1000-500 mb thickness prognoses to obtain a 
1000-mb forecast. The U. S. Naval Fleet Numerical Weather Facility 
(FNWF) has utilized an empirically adjusted advective technique for 
thickness prediction. A more complete thickness forecast model includ- 
ing vertical motion, heating and friction is presently under investiga- 
tion by Professor G. J. Haltiner of the U. S. Naval Postgraduate School. 
This program utilizes a thickness tendency derived from the first law of 
thermodynamics and is shown in a subsequent section of this paper. In 
some of the earlier applications of this equation some of the terms have 
been omitted by way of simplification because of the lack of information 
necessary to calculate these terms, for example, diabatic heating. Re- 
cently FNWF has made computations of the heat exchange on a regular 
basis. With this data now available it is appropriate to consider 
whether a better estimate of a thickness tendency might be obtained if 


more terms of the equation were retained and fewer simplifying assumptions 


made. 

The purpose of this investigation is to conduct a first step in a 
determination of the usefulness of this equation, that is, to statisti- 
cally correlate the observed 1000-500 mb thickness changes with estimates 
of the various terms in the thermodynamic equation. As shown by Panofsky 
and Brier [8] linear correlation coefficients give a measure of associa- 
tion between two variables which does not depend upon any arbitrary choice 
of units by which the original variables were measured. From the correla- 
tion coefficients obtained here each term will be examined as to its re- 
lation to the observed thickness change and its significance as a pre- 
dictor. From these results the value of this equation in a prediction 


scheme can be estimated. 
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Hae AREA OF INVESTIGATION 

FNWF routinely analyzes data over the Northern Hemisphere on a 
63 x 63 grid, with a grid distance of 381 KM. 

As shown later in this paper the scheme used in obtaining upwind 
points for data extraction was designed for use on a 22 x 22 grid extract- 
ed from any section of FNWF's 63 x 63 grid. The total heat exchange 
calculations were originally derived for use over the oceans only, how- 
ever FNWF calculates this term over the whole grid, and then plots only 
the ocean points. Because the land area points were to be kept at a 
minimum, all available values of the heat exchange were used. The known 
decrease in the accuracy of the value over land as compared to that over 
the ocean nevertheless was accepted in this study and recognized as a 
possible source of error. 

Due to the facts mentioned above two different areas of 22 x 22 grid 
size were selected for this study. These areas were chosen because most 
of the grid is oceanic with a minimum of continental area. 

Area 1, which will be referred to as the Atlantic area, is shown in 
figure 1. This area includes the North Pole and the Gulf Stream and has 
a minimum of equatorial area; the latter is near the boundary o*% Fi 's 
grid where results are considered to be less accurate. 

Area 2, which will be referred to as the Pacific area, is shown in 
figure 2. This area is almost entirely oceanic, but does contain more 
of the equatorial or boundary area. These areas will be treated separ- 


ately and the correlation coefficients will be calculated independently. 


ua 


3 BACKGROUND 
The thickness prediction equation may be derived in the following 


way. The first law of thermodynamics can be written 


y = oats t ip ee 
= Solute HMI Stale tll alia x UW) 
sf ii 1ii iv V (1) 


Combining the equation of state and the hydrostatic equation gives 


 — dz 
os tf Tp (2) 


Substituting (2) into (1) in terms (ii) and (iii) gives 


dz 
J 
- =~ 444,007) of V-V(55) + ee ty (3) 


Rearranging, we have 


29) « ¥y-9 (38) + Ru (35 -EE)- 8, W) 


Integrating between 1000 and 500 mb 


se eof v(55) 4p +[ Twey - (#8 aft (5) 


1600 1000 


where 
c= Be GE- 8) ° 


Now integrating over a 12-hour time period gives 


Nes {(w = +ff “o™. im (7) 


Oo /000 /000 / jas 
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The purpose of this investigation is to correlate the 12-hour thick- 
ness changes with numerical approximations of the terms constituting the 
right hand side of equation 7. Thus Ah will be correlated with 
estimated mean values of the horizontal advection, A, the vertical motion 
term, W, and the total heat exchange, H, for the 12-hour period. This 


relationship may be written 


Ahk = A +WwW +H re 


The coefficients of the terms in equation 7 can be considered as 
merely constants of proportionality i thus can be omitted in the 
statistical evaluation. 

The parameter 0 as defined here differs slightly from the value 
of GO used by FNWF. However, here again the difference can be con- 
sidered essentially as a constant of proportionality and will not affect 


this study. 
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4, DATA COLLECTION 

The basic data for this study was obtained for the three-day period 
9 to 11 February 1966 from the Fleet Numerical Weather Facility. 

The initial fields utilized were the SD-500 field which is the dis- 
turbance component (perturbation) of the 500-mb height field, the SR-500 
field which is the residual component (basic flow) of the 500-mb height 
field, the SD-1000 field which is the disturbance component at the 1000 
mb level and the SR-1000 field which is the basic flow at the 1000-mb 
level. 

These fields were used in a routine originally designed for the 
interpretation of meteorological satellite observations as developed by 
Nagle, Clark and Holl [7] and modified by Clark for use in obtaining the 
advective change in thickness, A. The scheme entails generation of 
stream function fields, Y , figures 3 and 4, from the SD and SR fields 
averaged over the 1000-500 mb layer. By a timewise interpolation of the 
12-hour distributions, hourly distributions of V were obtained which 
would act as evolving frame- of- reference steering fields. From these 

‘YY fields a horizontal velocity distribution was computed. This wind 
was based on the geostophic approximation and in this case was the mean 
wind in the layer. 

At this time the 12-hourly winds for a 22 x 22 subsection of the 
63 x 63 grid were extracted. This subregion was then interpolated in 
both time and space to give hourly winds into an expanded 44 x 44 array 
based on the 22 x 22 section. This maneuver quadruples the number of 
data points contained in the region. 


By utilizing the stream function and the wind field, an upstream 
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trajectory was projected and the position of a column of air at time 


th - 12 was determined. This was done for each data point in the region. 
Checks in the program flagged any point whose upstream position was out 

of the 44 x 44 grid. This point was then discarded for correlation 
purposes. 

Several steps were involved in the procedure used to determine A. 
Thickness fjelds, figures 5 and 6, were obtained from FNWF's D-500 and 
D-1000 fields. These D fields give the difference between the height 
of the surface (500 and 1000 mb, respectively) and its corresponding 
standard height. The trajectory scheme was then applied to this field 
to yield a new parameter, h, which is the combined local change of thick- 
ness, Ah , and the estimated horizontal advection, A. Figures 7 and 
8 depict this parameter in the Atlantic and Pacific Py, 25s ee 
The strong gradient in the extreme western section seen in these and 
subsequent maps was due to the check-flag system discussed previously. 
These are the areas of inflow where the upstream position was off the grid 
causing a zero to be inserted. While the zero value was eliminated be- 
fore correlations were made, the point was plotted as such resulting in 
tight gradients where none actually exist. 

The local or observed 12-hour change in thickness, figures 9 and 
10, was produced by subtracting the field at time cs - 12 from the field 
at time t.. To get an estimate of A at each point, Ah was subtracted 
from h; this was not completed until the format of the data had been 
changed and the correlation program was in use. 

A further modification of the routine allowed W to be computed. 


This vertical motion term was obtained by taking the product of 9 and 


iS 


u) , the vertical motion in (x,y,p,t) coordinates, and averaging in 
both time and space. 

To accomplish this sigma-l and sigma-3, the static stability for the 
layers, 1000-750 and 600-750 mb respectively, wete utilized. These 0 's 
were derived from a model by Holl, Bibbo and Clark [2] which is charact- 
erized by pressure levels defining a number of layers within which the 
static stability is independent of pressure. 

These terms were multiplied by the W fields at 850 and 500 mb. 
FNWF computes these parameters by a numerical solution of a quasi- 
geostrophic omega equation, as shown by Haltiner, Clarke and Lawniczak 
[1]. The resulting fields were TW g6q) figures 11 and 12, and TW eos 
figures 13 and 14. The trajectory scheme was again used to obtain up- 
wind values of these terms at time t. - 12 which were then averaged with 


values at t, to yield the 12-hour averages, 74) and CW cog? shown 


850 


in figures 15, 16, 17 and 18. The final averaging of CW and OW) 


850 500 
to obtain W was completed by the correlation program itself. 

The same procedure was used to obtain H. This term was derived from 
the fields of Qi the total air-ocean interchange as produced by FNWF 
shown in figures 19 and 20. The computation program of the total heat 
exchange by numerical methods and its adaption to the computer evolved 


from earlier manual methods developed by Laevastu [4]. The heating func- 


tion was obtained from the following equation: 
Qn = Q.- @, - Qo - Gr tee (9) 


where all Q's will be expressed in units of g-cal aie (24 yap ae 
Because of the relative newness of the computerized form of this product 
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each term will be briefly discussed. 


The insolation, Qs, is calculated from the equation 
3 
emer oT. «=, | ao,ccocc (10) 


which takes into account the total time that the sun is above the horizon, 
tae the noon altitude of the sun, A and a cloudiness factor, C. Since 
it is virtually impossible to take into consideration various types of 
clouds and their different effects on the insolation, Laevastu [3] 
devised an average cloudiness factor which requires only a knowledge of 
total cloud tenths. The results give a reasonable estimate of the total 
incoming radiation. 

The albedo, Q. is an estimate of the radiation reflected from the 
surface. Laevastu [5] statistically determined an empirical formula for 


calculating this term. 


2 
Qp= O15 Q, 7 CO-01 Q,) (11) 


This formula is valid only for daily calculations. A simplified formula 


for the computation of short term (hourly or 3-hourly) albedo is 


Qe = Q. =22 ol = Average solar altitude (12) 


The effective back radiation, Q» is the long wave radiation given 


off by the sea surface. This represents a heat loss at the surface 
Qe G2 oitamalnSCoTveaiOn nee) ilaluciOnone Sues) (13) 


and makes use of a linear formula of Lonnquist [6] which is corrected by 


Moeller's cloudiness factor. 


Wee 


The convective transfer of sensible heat, Qa? was obtained from 


Q, = 3900.26 +07 V\(Ww ae) > aha (14) 


) 


when heat is transferred from the water to the air; and 
6 (Ak oas WC) aati maces (15) 


when heat is transferred from the air to the water. 

The latent heat transfer, Q.> is that heat addition or loss due to 
evaporation or condensation. Lavastu [5] shows that a modified formula 
of Rohwer has been found to be most accurate in estimating evaporation; 


namely, 
Qe * (0:26 +0.77V)( 0.98 Cw ~ Ca), if O.98C,>@, (16) 


A modified form of this equation can be used for condensation at the sea 


surface, 
Qe = 9077 V (0:98 Gy — Ca) Le , if 0.9 <& (17) 


To calculate these terms and get a total heat exchange product 
certain basic synoptic numerical analyses must be available. These are: 
Surface winds, surface air temperature and water vapor pressure and 
clouds. These are all routinely produced at FNWF. 

By utilizing the trajectory program again, the to - 12 values of 
Q were obtained. Averaging over the 12 hours produced H, figures 21] 


and 22. No further data was needed. 
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op OBJECT AND METHOD INVESTIGATION 

The purpose of this investigation was two-fold. An attempt was 
made a) to statistically correlate the observed change in thickness 
with estimated values of the terms on the right hand side of equation 
(7) and b) to estimate the significance of each term in a prediction 
scheme. 

To accomplish this a multiple regression and correlation analysis, 
the BMDO3R programmed for the CDC 1604 electronic computer and converted 
for COOP MONITER Fortran 63, was utilized. This program will perform 
linear-multiple correlation analysis on data in any combination by allow- 
ing the designation of a dependent variable and deletion of sets of in- 
dependent variables. This can be done as many times as desired. 

At the completion of the data collection phase there were five 
parameters stored for each data point over a total of 8229 Atlantic and 
6739 Pacific points which were considered good for correlation purposes. 


These variables were: (1) hr s (2) TW (3) TW 43 (4) H ; and OyAh. 


850° 
To place this data in a form appropriate for the BMD program the format 
had to be changed from a row-wise listing to a column presentation. 
After accomplishing this, spot checks were made to ensure that the data 
had not gotten out of phase during the transformation. 

Since the parameters A and W had not actually been obtained as yet, 
the transgeneration feature of BMDO3R was used. This allows for the 
generation of new variables by performing certain mathematical opera- 
tions on the original variables. Here A, variable 6, was obtained by 
the subtraction (1) - (5) = (6) and W, variable 7, was produced from 


[(2) + (3)]% = (7). The data was then ready for the correlation analy- 


sis. 
19 


6. RESULTS AND DISCUSSION 

The data was divided into Atlantic and Pacific regions with identi- 
cal calculations made for each area. Initially, correlation analyses 
were carried out to show relationships between pairs of variates. Many 
combinations of independent and dependent variables were used. The re- 
sults are shown in Tables 1 and 2. 

After examining the figures, multiple linear correlations were made. 
This was done by progressively adding the term from equation (8) with the 
highest partial correlation coefficient. 

As a side operation, a multiple correlation analysis was made on 


the terms of the equation 


Rc Wt F 
These results are shown in Tables 3 and 4. 

The values of R obtained from the Atlantic data were slightly higher 
than those in the Pacific in almost every case. This can be partially 
explained by a) the scarcity of reporting stations in the Pacific result- 
ing in larger interpolated areas, and b) the Pacific area contains a 
larger band of equatorial or boundary region where the analysis is known 
to be less accurate. 

The coefficient obtained in the linear correlation of Ah and H 4% 
“extremely low in both the Atlantic and Pacific. 

The coefficient resulting from the correlation of Ah and W was lower 
than expected but was approximately the same in both Atlantic and Pacific. 

Ahk was correlated with TW eo and TW) 549 separately to see if it was 


necessary to obtain W or if the value of the parameter obtained at a 


single level would be sufficient. In the Atlantic there is little 
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Dependent Independent 
Variable Variable 
Ah A 
Ah Ww 

Ah. : 
4D he TWaso 
Ah TWeoo 
Ww A 
Ww -) 
A H 
h W 
ae Ww 


0099 


. 1936 


.0655 


.1729 


. 1850 


»4513 


. 2200 


a7 99 


. 2425 


. 3802 


7394 


.0375 


. 0043 


.0299 


0342 


. 2037 


0484 


.0324 


.0588 


. 3366 


Table 1. Results of non-multiple linear correlation analyses, 
Atlantic. R is the correlation coefficient and R2 is the co- 


efficient of determination. 
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Dependent 
Variable 


TAN oe 


Ah 


| 


A > 


°° 


Table 2. 


Independent 
Variable 


A 


= 


.7967 


e752 


.0985 


. 2440 


.0825 


. 2569 


0404 


.0962 


.0483 


BOIL 


6347 


.0302 


.0097 


0598 


.0068 


.0660 


.0016 


.0093 


.0023 


. 2886 


Results of non-multiple linear correlation analyses, 


Pacific. 
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Dependent Independent R R 
Variable Variable 
Ah YX 8599 7394 
Ah A +tWwW .8871 . 7869 
Ah A +W+H 8889 7901 
KA w+ 5920 3505 
Table 3. Results of the multiple-linear correlation analysis, 
Atlantic. © 
Dependent Independent R R? 
Variable Variable 
N\A A .7967 6347 
Ahk A +W 8562 7331 
Ah Atweit 8625 7439 
oO 
h Wt 5418 2936 


Table 4. Results of the multiple-linear correlation analysis, 
Pacific 
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difference in the coefficients but in the Pacific, the three differ 
considerably. 

The multiple correlation coefficient in the Atlantic increased 
from .8599 to .8871 to .8889 as W and H were added to A. The increase 
in the Pacific was somewhat greater, going from .7967 to .8562 to .8625 
as these terms were added. Here the inclusion of the vertical motion is 
more significant. 

The analysis of with W and H was attempted with the thought 
that hh might be used in a cloud forecast scheme currently under study. 
Although the correlations in both areas are significant, they are certain- 
ly not high enough to consider replacing vertical velocities by A in a 
cloud forecasting scheme. 4 


The remaining analyses clearly show no significance and will not be 


discussed further. 
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7 CONCLUSIONS AND RECOMMENDATIONS 

As a result of this study, which was made using wintertime data 
only, the following conclusions were reached. 

(1) Observed changes in thickness over a 12-hour period are due 
almost entirely to horizontal advection, at least as approximated in this 
study. This result lends strong support to the simple advective scheme 
for predicting thickness which has been used operationally by FNWF dur- 
ing the past several years, i.e. 

DN pee, eee acable Ns 

(2) The vertical motion term, W , though of secondary importance, 
does contribute significantly to the change in thickness and should be 
retained in a prediction scheme for increased accuracy. Of course, a 
decision must be made as to whether the improvement in forecasting skill 
is worth the time and effort necessary to include vertical motion. 

(3) The heat exchange term, H, as presently being calculated at 
FNWF and as utilized in this study contributes virtually nothing to the 
change in thickness and should be omitted in a thickness prediction scheme. 

As pointed out earlier, this study dealt only with winter ‘data. An 
obvious and necessary extension should be the application to data from 
the other seasons, especially summer. 

It is hoped that the results of this study will provide guidance to 
individuals or groups engaged in the development of ceWethess prediction 


schemes. 
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Figure 20 ‘ Total heat exchange, 


68 





a { 
¢ 


m a z : a“ 3 
2 : . 


SE 04 1 ; z e e 7 a a 2 = 13 . “ . ~ 
~ t= 








2tyTOeg ‘H ‘w190q SulRZeoH ueoW 
» 
KLE e nT a at 


es ¥ e 

* 

. * 

. 
° ° 
a . 
* 
° 
e ° 


ZZ eAnsTy ' rs 


4 : a 

a 
al 
ioog: 
Pipe. 


®e 
ate C- Ds Be fe, 
. cits ar fe by 
° io ~ . ., ° a 
ae oe ar Br ass 
es 2 *. e 
. ee ee bien : ° 
.. ce. 
ae Se Bi 
ae oy 
tee _* 
OE tus 





10, 


Lae 


INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library 
U. S. Naval Postgraduate Schoo 
Monterey, California 


LIT Wayne R. Lambertson 
3906 Ocean Blvd. 
Brant Beach, New Jersey 


Department of Meteorology & Oc 


U. S. Naval Postgraduate School 


Monterey, California 


Office of the U. S. Naval Weat 
U. S. Naval Station (Washingto 
Washington, D. C. 20390 


Officer in Charge 

Naval Weather Research Facilit 
U. S. Naval Air Station, Bldg. 
Norfolk, Virginia 23511 


Commanding Officer 

FWC / JTWC 

COMNAVMAR 

FPO San Francisco, Calif. 966 


Commanding Officer 
U. S. Fleet Weather Central 
FPO Seattle, Washington 98790 


Commanding Officer 
U. S. Fleet Weather Central 
FPO San Francisco, California 


Commanding Officer 
U. S. Fleet Weather Central 
FPO New York, New York 09540 


Commanding Officer 
Fleet Weather Central 
Navy Department Washington, D. 


No. Copies 
20 
7 
1 
1 
eanography 1 
her Service l 
mn Navy Yard Annex) 
1 
y 
R-48 
1 
30 
1 
1 
96610 
1 
1 


C. 20390 


1% 


eZ. 


i. 


14, 


ES 


16. 


iy. 


ie . 


i. 


20. 


ZA. 


22. 


No. Copies 


Commanding Officer 

Fleet Weather Central 

U. S. Naval Air Station 
Alameda, California 94501 


Commanding Officer and Director 
Navy Electronics Laboratory 
Attn: Code 2230 

San Diego, California 92152 


Officer in Charge 

Fleet Numerical Weather Facility 
U. S. Naval Postgraduate School 
Monterey, California 


U. S. Naval War College 
Newport, Rhode Island 02844 


Director, Naval Research Laboratory 
Attn: Tech. Services Info. Officer 
Washington, D. C. 20390 


Office of Chief Signal Officer 
Research and Development Division 
Department of the Army 
Washington, D. C. 


Commander 

Air Force Cambridge Research Center 
Attn: CROOTR 

Bedford, Massachusetts 


Geophysics Research Directorate 
Air Force Cambridge Research Center 
Cambridge, Massachusetts 


Program Director for Meteorology 
National Science Foundation 
Washington, D. C. 


American Meteorological Society 
45 Beacon Street 
Boston, Massachusetts 


Commander, Air Weather Service 
Military Airlift Command 

U. S. Air Force 

Scott Air Force Base, Illinois 62226 


73 


— a 
a ° 


Zon 


24. 


25: 


26. 


27. 


26" 


29% 


30. 


be 


O7s 


33. 


34. 


a 


No. 


U. S. Department of Commerce 
Weather Bureau 
Washington, D. ©. 


Office of Naval Research 
Department of the Navy 
Washington, D. C. 20360 


Department of Meteorology 
University of California 
Los Angeles, California 


Department of the Geophysical Sciences 
University of Chicago 
Chicago, Illinois 


Department of Atmospheric Science 
Colorado State University 
Fort Collins, Colorado 


Department of Engineering Mechanics 
University of Michigan 
Ann Arbor, Michigan 


School of Physics 
University of Minnesota 
Minneapolis, Minnesota 


Department of Meteorology 
University of Utah 
Salt Lake City, Utah 


National Center for Atmospheric Research 
Boulder, Colorado 


Department of Meteorology & Climatology 
University of Washington 
Seattle, Washington 98105 


Department of Meteorology 
University of Wisconsin 
Madison, Wisconsin 


Department of Meteorology 
Florida State University 
Tallahassee, Florida 


Department of Meteorology 
Massachusetts Institute of Technology 
Cambridge, Mass. 02139 


74 


Copies 


2 


36. 


37. 


38. 


Bole 


40. 


41. 


42. 


43. 


44. 


45. 


46. 


No. 


Department of Meteorology 
Pennsylvania State University 
University Park, Pennsylvania 


Hawaii Institute of Geophysics 
University of Hawaii 
Honolulu, Hawaii 


University of Oklahoma 
Research Institute 
Norman, Oklahoma 


Atmospheric Science Branch 
Science Research Institute 
Oregon State College 
Corvallis, Oregon 


The University of Texas 

Electrical Engineering Research Laboratory 
Engineering Science Bldg. 631A 

University Station 

Austin, Texas 78/712 


Department of Meteorology 
Texas A&M University 
College Station, Texas 77843 


Lamont Geological Observatory 
Columbia University 
Palisades, New York 


Division of Engineering and Applied Physics 
Room 206, Pierce Hall 

Harvard University 

Cambridge, Massachusetts 


Department of Mechanics 
The Johns Hopkins University 
Baltimore, Maryland 


University of California 
E. O. Lawrence Radiation Laboratory 
Livermore, California 


Department of Astrophysics & Atmospheric Physics 


University of Colorado 
Boulder, Colorado 


75 


Copies 


1 


47. 


48. 


49. 


50. 


Sle 


Dee 


53) 


54. 


ae 


56. 


Weather Dynamics Group 
Aerophysics Laboratory 
Stanford Research Institute 
Menlo Park, California 


Meteorology International, Inc. 
P.O. Box 1364 
Monterey, California 93940 


The Travelers Research Center, Inc. 
650 Main Street 
Hartford, Connecticut 


Meteorology Section 

Aracon Geophysics Co. 
Virginia Road 

Concord, Massachusetts 01742 


Librarian 

GCA Technology Division 

GCA Corporation 

Bedford, Massachusetts 01730 


Department of Meteorology 
University of Melbourne 
Grattan Street 

Parkville, Victoria 
Australia 


Bureau of Meteorology 
Department of the Interior 
Victoria and Drummond Streets 
Carlton, Victoria 

Australia 


International Antarctic Analysis Centre 
468 Lonsdale Street 

Melbourne, Victoria 

Australia 


C.S.1.R.0. 

Division of Meteorological Physics 
Station Street 

Aspendale, Victoria 

Australia 


Institut f£. Meteor. u. 

Geophysik Universitat Innsbruck 
Schopfstrassee 41, Insbruck 
Austria 


76 


No. Copies 


i 


oe 


58. 


a9 . 


60. 


61. 


62. 


63. 


64. 


65. 


66. 


67. 


Department of Meteorology 
McGill University 
Montreal, Canada 


Central Analysis Office 
Meteorological Branch 
Regional Adm. Building 
Inter. Airport 

Dorval, Quebec, Canada 


Meteorological Office 
315 Bloor Street West 
Toronto 5, Ontario, Canada 


Department of Meteorology 
University of Copenhagen 
Copenhagen, Denmark 


Institute of Meteorology 
University of Helsinki 
Helsinki - Porthania, Finland 


Institut fur Theoretische Meteorologie 
Freie Universitat Berlin 
Berlin-Dahlem 

Thiel-allee 49 

Federal Republic of Germany 


Meteorological Institute 
University of Thessaloniki 
Thessaloniki, Greece 


Meteorological Service 
44, Upper O'Connell Street 
Dublin 1, Ireland 


Department of Meteorology 
The Hebrew University 
Jerusalem, Isreal 


Geophysical Institute 
Tokyo University 
Bunkyo-ku 

Tokyo, Japan 


Meteorological Research Institute 


Kyoto University 
Kyoto, Japan 


77 


. Copies 


68. 


69. 


70. 


(Oke 


123 


73. 


74, 


Pe 


76. 


ee 


78. 


No. 


Department of Astronomy & Meteorology 
College of Liberal Arts and Sciences 
Seoul National University 

Tbng Soong Dong, Chong No Ku 

Seoul, Korea 


Central Meteorological Office 
I Song Wul Dong, Sudaemon Ku 
Seoul, Korea 


Department of Meteorology 
Instituto de Geofisica 
Universidad Nacional de Mexico 
Mexico 20, D. F., Mexico 


New Zealand Meteorological Service 
P.O. Box 722 
Wellington, G. E. New Zealand 


Institutt for Teoretisk Meteorologi 
University of Oslo 
Blindern, Oslo, Norway 


Institute of Geophysics 
University of Bergen 
Bergen, Norway 


Pakistan Meteorological Department 
Insittute of Meteorology & Geophysics 
Karachi, Pakistan 


Royal Swedish Air Force 
M.V.C. 
Stockholm 80, Sweden 


Department of Meteorology 
Imperial College of Science 
South Kensington 

London S. W. 7, United Kingdom 


Meteorological Office 
London Rd. 

Bracknell 

Berkshire, United Kingdom 


National Research Institute for Mathematical Sciences 
Ges. Te. 

P.O. Box 395 

Pretoria, Union of South Africa 


78 


Copies 


79: 


Mr. James R. Clark 

Meteorology International, Inc., 
419 Webster St., 

Monterey, Calif. 93940 


79 


No. Copies 


I 





Unclassified 


Security Classification 


DOCUMENT CONTROL DATA - R&D 


(Security claseitication of title, body of abstract and indexing annotation muet be entered when the overall repart te Claeailied) 
1. ORIGINATING ACTIVITY (Corporate author) 24. REPORT SECURITY CLASSIFICA TION 
U. S. Naval Postgraduate School Unclassified 
Monterey, California 


3. REPORT TITLE 


A Correlation Analysis of the 1000-500 mb Thickness Tendency 


4. DESCRIPTIVE NOTES (Type of report and inclueive detea) 
None 


8. AUTHOR(S) (Last name, rat name, initial) 





LAMBERTSON, Wayne R., LT, USN 


6. REPORT DATE Ja. FOTAL NG. OF PAGES 7b. NO. OF ARFS | 
June 1966 80 9 


Ga. CONTRACT OR GRANT NO. 9a. ORIGINATOR’S REPORT NUMBER(S) 


N.A. 





b. PROJECT NO. N A 


9b. OTH PORT 3) (A 
ah ZR oh ° NO(S) (Any other numbere thet may be ecoigred 


N.A. 
10. AVAIL ABILITY/LIMITA TION NOTICES 


11. SUPPL EMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 
. Chief of Naval Operations (OP-09B7) 
None Navy Department 
Washington, D. C. 20350 
- ABSTRACT 
A correlation analysis of a thickness tendency equation including heating 
and vertical motion terms was carried out. One of the main purposes of the 


investigation is to provide a firmer foundation for a 1000-500 mb thickness 


forecast scheme. 


A large data sample was used and single and multiple linear coefficients 
were calculated. These were computed for a wintertime situation with separate 
analyses for Atlantic and Pacific areas. 

The resulting correlation coefficients for the advection, vertical motion 
and heating terms gave substance to the hypothesis that the observed change in 
thickness may be fairly well approximated by the change due to horizontal 
advection only. The vertical motion term gave a smaller but significant corre- 
lation, while heat exchange terms did not correlate significa... .y. 





DD ae 1473 Unclassified 


81 Security Classification 


Unclassified 


~~ Security Classification 





Correlation Analysis 

1000-500 mb thickness tendency 
thickness 

thickness forecast scheme 

heat exchange 


LINK A LINK B LINK C 
KEY WORDS 


INSTRUCTIONS 


1, ORIGINATING ACTIVITY: Enter the name and address 
of the contractor, subcontractor, grantee, Department of De- 
fense activity or other organization (corporate author) issuing 
the report. 


2a. REPORT SECURITY CLASSIFICATION: Enter the over 
all security clasaification of the report. Indicate whether 
‘*Restricted Data’’ ia included. Marking is to be in accord 
ance with appropriate security regulations. 


26. GROUP: s.specified. » 
rective 5200. ite Coen yneninminers tr eae ora - 


Enter 
the group number. Also, when applicable, show that optional 
markings have been used for Group 3 and Group 4‘as author- 


ized. 


3. REPORT TITLE: Enter the complete report title in all 
capital letters. Titles in all cases should be unclassified. 
If a meaningful title cannot be selected without classifica- 
tion, show title classification in all capitals in perenthesis 
immediately following the title. 


4. DESCRIPTIVE NOTES: If appropriate, enter the type of 
report, e.g., interim, progress, Summary, annual, or final. 
Give the inclusive dates when a specific reporting period is 
covered. 


5. AUTHOR(S): Enter the name(s) of author(s) as shown on 
or in the report. Enter last name, first name, middle initial. 
If military, show rank and branch of service. The name of 
the principal author is an absolute minimum requirement, 


6. REPORT DATE: Enter the date of the report as day, 
month, year; or month, year. If more than one date appears 
on the report, use date of publication. 


7a. TOTAL NUMBER OF PAGES: The total page count 
should follow normal pagination procedures, Le., enter the 
number of pages containing information 


7b. NUMBER OF REFERENCES. Enter the total number of 
references cited in the report. 


8a. CONTRACT OR GRANT NUMBER: If appropriate, enter 
the applicable number of the contract or grant under which 
the report was written. 


8b, &, & 8d. PROJECT NUMBER: Enter the appropriate 
military department identification, such as project number, 
subproject number, system numbers, task number, etc. 


9a. ORIGINATOR’S REPORT NUMBER(S): Enter the offi- 
cial report number by which the document will be identified 
and controlled by the originating activity. This number must 
be unique to this report. 


9b. OTHER REPORT NUMBER(S): If the report has been 
assigned any other report numbers (either by the originator 
or by the sponsor), also enter this number(s). 


10. AVAILABILITY/LIMITATION NOTICES: Enter any lim- 
itations on further dissemination of the report, other than those 


DD 1473 (BACK) 


FORM 
1 JAN 64 


82 





aS by security classification, using standard statements 
such as: 





(1) ‘Qualified requesters may obtain copies of this 
report from DDC.’’ 
(2) ‘“‘Foreign announcement and dissemination of this. 
report by DDC is not authorized.’’ 
(3) “‘U. S. Government agencies may obtain copies of 
this report directly from DDC. Other qualified DDC 
ee users s hall requéatathrough jh)? ideas (8 eee 
99 
(4) ‘*U. S. military agencies may obtain copies of this . 
report directly from DDC Other qualified users 
shall request through 
iB 
(S) ‘*All diatribution of this report is controlled. Qual- 


ified DDC uaers ahalil request through 


8 
» 


If the report haa been furnished to the Office of Technical . 
Services, Department of Commerce, for sale to the public, indi- 
cate this fact and enter the price, if known 


11. SUPPLEMENTARY NOTES: Use for additional explana- 
tory notes. 


12. SPONSORING MILITARY ACTIVITY: Enter the name of 
the departmental project office or laboratory sponsoring (pay 
ing for) the research and development. Include address. 


13. ABSTRACT: Enter an abstract giving a brief and factual 
aummary of the document indicative of the report, even though 
it may alao appear elsewhere in the body of the technical re- 
port. If additional space ia required, a continuation sheet shall 
be attached. 


It is highly desirable that the abstract of classified reporta 
be unclassified. Each paragraph of the abstract shall end with 
an indication of the military security classification of the in- 
formation in the paragraph, represented as (TS), (S), (C), or (U). 


There is no limitation on the length of the abstract. How- 
ever, the suggested length is from 150 to 225 words. 


14. KEY WORDS: Key words are technically meaningful terma 
or short phrases that chsracterize a report and may be used as 
index entries for cataloging the report. Key words must be 
selected so that no security classification is required. Identi- 
fiers, such as equipment model designation, trade name, military 
project code name, geographic location, may be used as key 
words but will be followed by an indication of technical con- 
text. The assignment of links, rales, and weights is optional. 





Unclassified 
Security Classification 





7 
i 
, 
a 
- 
ah 
_ 7 7 
i 
= — — ae 
; a) 





a. & 











2) 

4 
2: = 
-~ —_— 


